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ABSTRACT 

Many observed massive star-forming z sa 2 galaxies are large disks that exhibit irregular morphologies, with 
sa 1 kpc, sa 10 Mg clumps. We present the largest sample to date of high-resolution cosmological SPH 
simulations that zoom-in on the formation of individual sa 10 10 5 M Q galaxies in sa 10 12 M Q halos at z sa 2. 
Our code includes strong stellar feedback parameterized as momentum-driven galactic winds. This model re- 
produces many characteristic features of this observed class of galaxies, such as their clumpy morphologies, 
smooth and monotonic velocity gradients, high gas fractions (/ g sa 50%) and high specific star-formation rates 
(> 1 Gyr -1 ). In accord with recent models, giant clumps (M c i ump «(5 x 10 8 - 10 9 ) Mg) form in-situ via gravi- 
tational instabilities. However, the galactic winds are critical for their subsequent evolution. The giant clumps 
we obtain are short-lived and are disrupted by wind-driven mass loss. They do not virialise or migrate to the 
galaxy centers as suggested in recent work neglecting strong winds. By phenomenologically implementing the 
winds that are observed from high-redshift galaxies and in particular from individual clumps, our simulations 
reproduce well new observational constraints on clump kinematics and clump ages. In particular, the obser- 
vation that older clumps appear closer to their galaxy centers is reproduced in our simulations, as a result of 
inside-out formation of the disks rather than inward clump migration. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: structure 



1. INTRODUCTION 

Star-forming disk galaxies at z sa 2 differ from their lo- 
cal counterparts in several aspects: they are more gas-rich 
and rapidly star-forming, and have irregular morphologies 
and high gas velocity dispersions. Spatially resolved obser- 
vations of their rest-frame UV and optical continuum light 
often reveal large (sa 1 kpc) clumps (lElmegreen et al.l [20051 
l2009t IForster Schreiber et alJl201 lbl) . Clumpy morphologies 
are also observed wit h star-formation tracers such as Ha 
(e.g. the SINS survey; iGenzel et al]|2008l 1201 lb . The spec- 
trally resolved line emissions reveal ordered rotation and high 
velocity dispersions, widespread across the disks at sa 50 - 
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Theoretically, the large z sa 2 star-forming galaxies have 
been considered as marginally-unstable thick disks, where the 
large random motions set the large masses an d sizes of un- 
stable regions that collapse into giant clumps (llmmeli et alJ 
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iBurkert et alfeOlOtlCeverino et al.ll2010l) . In this picture, the 
high star-formation rates res ult from the high growth rates 
of cosmic structures at z sa 2 (iGenel et al.ll2008f) . as gas fun- 
nels directly from the cosmic web to the vicinity of galax- 
ies in a 'cold-mode' accre tion ( iKeres et al.l l2005; Ker es et alJ 
120091; iDekel et all l2009al) . However, several observational 
clues have not yet been addressed in this framework. First, 
vigorous outflows are observed from a vari ety of star-forming 
galaxies at high red s hift, both globally (e . g.lPettini et al.l2000fc 
IWeiner et all l2009t IShapiro et all 120091; ISteidel et alJ l2010h . 
and from individual giant clumps (IGenzel et al.ll201 II) . Sec- 
ond, high signal-to-noise observations of individual giant 
clumps reveal only minor kinematical signatures, sometimes 
indicating dynamical masse s lower than altern ative indepen- 
dent clump mass estimates dGenzel et alj|201 11) . Third, com- 
parisons of abundance and clustering of galaxies and dark 
matter halos suggests that at z sa 2, much like at z = 0, typ- 
ically not more than sa 10% of the bary ons associated with 
halos of any mass have turned into sta rs ( iMoster et alj|2010l ; 
iBehroozi et al.l2010tlWake et alj201 lb . causing some tension 
with the high efficiency of 'cold-mode' accretion. 

Further on, one aspect of the theoretical picture that de- 
veloped in recent years is that the giant clumps lose angu- 
lar momentum and energy due to dynamical friction in the 
disk, clump-clump interactions, and internal torques within 
the disk, and thereby spiral in to the center of the galaxy. This 
has been based on two main lines of argument. First, age 
gradients have been inferred in clumpy disks from the pres- 
ence of a redder/o lder component often coinciding or close to 
the ga laxy center ( Ivan den Bergh et alj!996t lElmegreen et alJ 
120091) . Also, a correlation has been found between the rel- 
ative dynamical dominance of that central component to a 
(metallicity-based) age estimate for the system, which has 
been interpreted as indicative of a time sequence, where older 
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systems have had more time to tr ansform part of the ir disks 
into bulges via clump migration dGenzel et al.ll2008l) . Sec- 
ond, numerical simulations, both of isolated unstable disks 
dlmmeli et alj|2004at iBournaud etaT] 120071: lElmegreen et alj 
l2008t i Bournaud et al. 2009) a nd of cosm ological z ~ 2 disks 
(lAgertz et alJl2009t ICeverino et al.ll2010h . find the clumps to 
evolve in that fashion, and theoretical considerations also 
sugge st there should be a significant inflow rat e inside the 
disk (iDekel et al.1 l2009bt IBournaud et all 1201 ll) . Recently 
iForster Schreiber et al.l (1201 lbl) used deep higher-resolution 
NIC2 imaging to identify more robustly clumps in rest-frame 
optical emission, and, in combination with AO-assisted SIN- 
FONI Ha observations and ACS imaging, radial trends in 
clump stellar age with distance to galaxy center were inferred. 
This result appears consistent with the the oretical inspiraling 
picture, as noted by lCeverino et al.l (1201 ll) . 

Therefore, a crucial issue is whether giant clumps can sur- 
vive long enough to become stellar-dominated and/or sink 
to the centers of their galaxies, contributing to the buildup 
of bulges. Conflicting theoretical claims hav e been made 
dMurrav et al.l 120101: iKrumholz & Dekell2010l) , while prob- 
ing that question directly with observations is a difficult task. 
Therefore, the qualitative agreement between simulations and 
observations has been interpreted as supporting the picture 
of long-lived migrating clumps. However, none of the cos- 
mological or iso lated simulations reported so far (but see 
ISales et aTll2010l) . which find long-lived migrating clumps, 
include strong galactic winds, which as noted already are di- 
rectly observed. As the 'clumpy phase' appears to be ubiq- 
uitous at z ps 1 — 3, it is of fundamental importance for our 
understanding of galaxy formation to know what role giant 
clumps play in building bulges/spheroids. 

In this paper, we demonstrate as a proof of principle that 
incorporating galactic superwinds resulting from stellar feed- 
back into the existing theoretical picture, can have important 
implications for our understanding of z ~ 2 disks and their 
evolution. We use hydrodynamical cosmological simulations 
to investigate the formation and properties of a sample of star- 
forming disks. We use a phenomenological model for the gen- 
eration of galactic winds to test a scenario in which galactic 
superwinds are a major driving force in shaping galaxy disks 
at high redshift. The phenomenological a priori parameters 
of our model match the observed properties of winds in the 
type of galaxies we are considering. Our results are heav- 
ily dependent on this phenomenology, hence the uncertainty 
in the observations of such winds is a source of uncertainty 
for the basic assumptions of our model. However, we find 
that the implications of implementing such winds in simula- 
tions, namely that the giant clumps are short-lived transient 
features, and thus do not survive to migrate to the centers of 
their galaxies to form a bulge, are consistent with the existing 
observational constraints. Hence, we propose a new scenario 
for the evolution of z ~ 2 clumpy disks. 

2. THE SIMULATIONS: CODE AND SETUP 

We run 'zoom- in' cosmological simulations focused on in- 
dividual halos taken from the 72/T 1 Mpc cosmo logical dark 
matter simulation presented in lOser et all (1201 Oh for the fol- 
lowing cosmological parameters: Q, m = 0.26, J7a = 0.74, f2j = 
0.044, h = 0.72, n = 0.95 and <j 8 = 0.77. As we focus on 
massive disks with high star-formation rates (SFR), we se- 
lect at z = 2 halos of ps 10 12 M Q with instantaneous dark mat- 
ter growth rates exceeding 500 M© yr" 1 and that had no ma- 
jor merger (mass ratio < 3 : 1) since z = 3 (see iGenel et ail 



120081) . The second criterion selects approximately half of the 
halos of that mass at that redshift, and together these crite- 
ria are met by ps 15% of those halos, i.e. our galaxies form 
in halos that are neither the most 'typical' nor very 'special' 
or rare. From this sample, we randomly selected nine halos 
and generated 'zoom-in' initial conditions for r e-simulations 
including baryons, as in lOser etal.ll20Toll2011l Our scheme 
results in high-resolution regions of ps 5 Mpc comoving, on 
a side, i.e. approximately 1000 times larger (in volume) than 
our main halos at z ps 2. The mass resolution of our simula- 
tions is 8 x 10 5 M Q and 5 x 10 6 M Q for baryonic and high- 
resolution dark matter particles, respectively. The gravita- 
tional softening lengths of baryonic particles are 200/j -1 pc, 
constant in comoving units, resulting in physical softening 
lengths of pa 90 pc at z = 2. 

The cosmological box is evolv ed to z = 2 using a version of 
the N-body/SPH code Gadset- 2 dSpringelll2005l) that i s very 
simila r to the one described in Oppenh eimer & Pavel d2006l 
12008b . This version includ es ionisation and heating by a uni- 
form background radiation dHaardt & Madaul2001l) in the op- 
tically thin limit, atomic cooling down to T = 10 4 /T from hy- 
drogen and metals, star-formation and feedback, as well as 
mass loss and metal enrichment from AGB stars and super- 
novae of types II and la. 

The lOppenheimer & Pavel d2008l) version of the code uses 
the ISpringel & Hernquisll (l2003h star-formation prescription, 
where star s form stochastically from the gas accor d ing to 
a ISchmidtl d!959l) law that reproduces the iKennicuttl dl998l) 
gas surface density to star-formation relation. In addition, 
that prescription includes a sub-grid model for supernova 
feedback that results in an effective equation of state for 
the star-forming gas P oc p 7(p) , where P is the gas pres- 
sur e, and p its physical density . The running index j(p) in 
the ISpringel & Hernquisll d2003l) model is rather high at low 
and medium densities (a stiff equation of state), but becomes 
negative at high densities, such that the thermal Jeans mass 
Mj can rapidly drop below the resolution of the simulation 
dHopkins & Ouataertll2010l) . The stiffness of this equation 
of state at medi um densities allows construction of stable 
disk models (e.g. ISpringel et al1l2005l) . It is less relevant for 
simulations of marginally-unstable di sks. Therefore, we re- 
placed the star-formatio n model in the lOppenheimer & Pavel 
(120081) code with the ISchaye & Palla Vecchial d2008h ap- 
proach. There, the usual cooling and heating operate at 
p < Peos, while for p > peos a polytropic equation of state 
with 7 = 4/3 is implemented. This ensures that the (ther- 
mal) Jeans mass Mj remains constant with varying density 
and does not become unresolved. Here peos = 0.1 cm -3 is de- 
fined such that Mj ps 1 .4 x 10 7 M Q at p = p E0S and T = 10 4 K. 
Star-formation is parameterised as 

^ = 7^)°' 26 , (1) 

dt lGyr p th 

with the star-formation density threshold p t h = p eos = 
0.1 cm -3 . This reproduces the | Kennicutt) (| 19981) star- 
formation relation (ISchaye & Palla Vecchiall2008l) . A further 
addition we introduced to the code is that we calculate the lo- 
cal mini mum (over all direct ions) of the gas surface density 
p 2 /\Vp\ dGnedin et al.ll2009h . and suppress star-formation if 
p 2 /\V p\ < 10 21 cm" 2 . This modification has only minor con- 
sequences for our results. 

The feedba ck scheme developed by lOppenheimer & Pavel 
d2006l 120081) builds on th e kinetic feedback scheme of 
Sprin gel & Hernquisll d2003l) . where star-forming gas parti- 
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cles are stochastically kicked with velocities in the v x a di- 
rection, where v is the particle's velocity prior to the ki ck and 
a its acceleration. lOppenheimer & Pavel d2006t 12008) intro- 
duced a significant modification to this mechanism by tuning 
the two parameters that control the wind. In their prescrip- 
tion, the magnitude of the kick is v w m = c(2 + 3\//l - 1) and 
the mass-loading factor 



where fi, jw 1 — 2 is the luminosity factor, a is the 've- 
locity dispersion' of the galaxy that is calculated from its 
mass 10 , and op is a consta nt (here ao = 300 kms" 1 , as in 
lOppenheimer & Pavel d2006b ). Introducing c-dependencies, 
these two parameters then scale with the mass of the galaxy, 
following the theory of momentum-driv en winds induced by 
the radiation pressure from young stars dMurrav et alj|2005t : 
IZhang & Thompsonl2010l) . 

The momentum budget availa ble from feedback from 
young stars is somewhat uncertain. iMurrav et alj d2005l) esti- 
mate the momentum injection from radiation pressure of OB 
stars (per unit stellar mass formed) to be w 200 kms -1 , and 
a similar amount of mom entum from supernova explosions. 
lOstriker & Shettvl d201 ll) . however, estimate the specific mo- 
mentum available from supernovae to be ps 3000 kms -1 . In 
addition, high dust column densities can provide large opti- 
cal depths that increase the amount of momen tum injection 
by ra diation pressure by factors of several (e.g. lMurrav et alj 
120101) . Finally, stellar winds can also inject similar amounts 
of momentum dLeitherer et al.lfl999l) . We use v W i„d = a(4 + 
4.29\//l— 1), such that the value that we use for the mo- 
mentum injection per unit stellar mass formed, r\\' Km d < 
2500 kms" , is a physically plausible value, even if on the 
high side. Note that we use higher wind velocities than in 
lOppenheimer & Pavel d2008l) . This ensures the escape of all 
wind particles from the disk, in particular from regions in the 
disk that have higher escape velocities than the mean one. 
This change is necessary, since we have a significantly higher 
resolution in this work, such that the center of the galaxy may 
become dense enough for its escape velocity to be higher than 
that at the edge of the disk. 

The scalings of the momentum-driven winds are in accor- 
dance with ob servational evidence for the scal ings of galac- 
tic winds (e.g. |MarH3l2Q^|^pkeietaU|2005l). The typical 
velocities we obtain with this model for the central galaxies 
considered are sa 400-700 kms -1 , in agreement with observa- 
tions. The typical mass-loading factor for our central galaxies 
is rj w 4, also in order-of-magnitude agreement with observa- 
tions, although the observational determination of the mass- 
loading fa ctors is far less certain than that of t he wind veloc- 
ities (e.g. lErb et alj|2006t iGenzel et al.ll2.01 ll and references 
therein). 

In the ISpringel & Hernquistl (120031) and 
lOppenheimer & Pavel d2008l) models, subsequent to be- 
ing given a velocity kick, gas particles are temporarily 
decoupled from hydrodynamic interactions to allow them 
to propagate out of their star-formation sites and eventually 

10 Galaxy dynamical masses are calculated based on an on-the- 
fly group finder, and then used to calculate a = 200[(M(jy n /5 ■ 
10 12 M0)4(i?(i)///o)]'' 3 . For calculation efficiency, a friends-of-friends 
(FOF) algorithm is used. We use a linking length b = 0.025(g^) '' /3 that 
is an improvement over the constant b used in Oppenheimer & Dave 1 2008), 
since it is calibrated to the SKID group finder at various redshifts. 



leave their galaxies. When simply giving gas particles 
momentum kicks rather than using a more sophisticated 
physical model (in a better-resolved interstellar medium), 
the only way to achieve winds that escape the vicinity of 
the galaxies (rather than efficiently becoming a galactic 
fountain) is by using the tem porary decoupling method 
dPalla Vecchia & Schavd 120081) . Moreover, even high 
resolution simulations of isolated galaxies have so far had 
difficulties in forming efficient winds when having full 
hydrodynamics and cooling (e.g. | MacLow & Ferraralll999l : 
iFuiita et all 120041: iPubois & Tevssierj|2008l) . However, such 
winds are observed, and eme rge out of z ~ 2 galaxies to large 
distances dSteidel et al.ll2010T) . Therefore, we keep the decou- 
pling prescriptions adopted by lOppenheimer & Pavel d2008l) . 
We have, however, modified the lOppenheimer & Pavd d2008) 
code to e nsure that the velocity kicks are perpendicular to the 
disk (e.g. lBordoloi et afll201 ll Bouche et al. in prep.). This is 
done by calculating the direction of the velocity kicks wind 
particles are given as (v— v ga i) x (a-a ga i), where v, a are the 
velocity and acceleration of the launched wind particle just 
before the kick, and v ga i, a ga i are the velocity and acceleration 
of the center-of-mass of the galaxy, i.e. v x a is calculated in 
the inertial frame of the galaxy. 

The theoretical considerations and actual observations of 
winds fro m galaxies are augmen t ed by further motivation for 
using the lOppenheimer & Pavel d2008l) model, namely this 
strong galactic winds model succesfully matches several im- 
portant observations: 

• Properties of the intergala ctic medium, such as rnetal- 
licity and ionization sta te d Oppenheimer & Davdl20M 
120091: Ibave et al.ll2?Tl0h . 

• The galaxy mass-metall icity relation dFinlator & Pavel 
l200a lDav6 et al.ll2011ab . 

• Hig h galactic g as fractions at high redshift (Section lXTI 
and iPave et alj|201 lab . 

• Low fraction of cosm ic baryons in galax ies at all red- 
shifts (Sectionl3~T1and lPave et aT]|201 lbl) . 

It should be kept in mind, however, that the simulations pre- 
sented in this work have much higher resolution than p revious 
works that employed the lOppenheimer & Pavel d2008h model. 
Therefore, some adverse side effects cannot be excluded be- 
fore a large cosmological box is run with a similarly high res- 
olution. Indeed, we corrected for one such side effect by using 
increased velocities (see above). 

3. RESULTS 

3.1. Global properties 

We begin by considering several global properties of the 
galaxies formed in our sample of nine 'zoom-in' simulations 
and comparing them with observational constraints. Table QJ 
gives halo masses, peak circular velocities, sizes, gas and stel- 
lar masses, gas fractions and SFRs for the main galaxies in 
our sample. This sample can be compared to a 'median disk' 
from the SINS survey, whose properties are given in Table [TJ 
as well. Bearing in mind the uncertainties in the observational 
quantities, which are of roughly a factor of two for the masses 
and SFRs, and the unavoidable differences in the definitions 
of some of the quantities between simulations and observa- 
tions (see footnotes to Table [TJ), it seems that our simulated 
galaxy sample is overall a good match to the disks observed 
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FIG. 1. — Global properties of nine 'zoom-in' simulations at z. = 2, where 
each color represents one simulation, with its main galaxy (circles) and the 
smaller galaxies that formed in neighbouring dark matter halos inside the 
zoomed-in region (dots). Only galaxies that are central in their own ha- 
los are shown, i.e. leaving out satellite galaxies that are already inside the 
main zoomed-in halo. In all panels s22 4 is marked with a filled circle. Top: 
SFR versus stellar mass. The observed Daddi et al. 1 2007) relation for z. = 2 
(solid; standard deviat ion is dashe d), Rodighiero et al. (2010) relation for 
1.5 < z < 2 (red) and IBauer et alj <2011l) relation for 1.5 < z < 3 (blue), 
all corrected for a[Chabrier 1 200^) initial mass function, are plotted on top 
of the simulation results, and demonstrate a reasonable match between our 
simulated galaxies and observations. Middle: Gas fractions (gas mass over 
baryonic mass) versus stellar mass. There is a trend of increasing gas frac- 
tions with decreasing mass, and the fractions for our main sample are f» 40%, 
in good agreement with observations at z, ~ 2.2 iTacconi et al. 2010, pluses). 
Bottom: 'Baryon conversion efficiency' as a function of halo mass. The 
highest efficiencies of es 10% are obtained at Mj m i > 1O 12 M0, in rough 
agreement with results obtained fr om matching observ ations to dark matter 
halo models, which are taken from Moster et al. 12010) f or z = 2.5 (dashed) , 
z = 1 .8 (dotted) and z = (dashed-dotted), as well as from lWake et alj <2011l) 
for 1 < z < 2 (solid). In all panels, the baryonic quantities (i.e. all but M200) 
are computed inside 3i?i/2, where R[/2 is the stellar half-mass radius within 
0.2R 20 o. 



in the SINS survey. Nevertheless, our simulated sample may 
be somewhat 'scaled-down' compared to the observed sam- 
ple, by a factor of approximately two in the masses and SFRs, 
and less than that in size. 

In Figure [T] panel (a) shows the SFR versus stellar mass 
(M») plane for central galaxies in the high-resolution region 
of each of our simulations. The simulated galaxies popu- 
late the sam e regio n as obse rved z ~ 2 galaxies, taken from 



Daddi et al.1 (120071) (bla ck), IBauer et al.l (1201 lb (blue) and 
Rodighiero et al.l (20101) (red). The scatter around the mean 
Dad di et al. (120071)" relation is roughly 0.5dex (dashed), and 



the uncertainty in those measurements is of similar magni- 
tude, as evident from the differences between the three lines. 
The SFR history in the central galaxies is rather constant, as 
demonstrated in Table Q]by comparing the current SFR to its 
past average. FigureQJb) presents gas fractions (/ g = M M ^ M ) 
versus stellar mass. There is a tre nd of increasing gas fractions 
with decreasing mass (see also iDave et all 1201 lal) . and the 
fractions for our main samp le are sa 40%, in go od agreement 
with observations at z ~ 2.2 (ITacconi et alj2010l pluses). Fig- 
ure [TJc) shows 'baryon conversion efficiencies' ( jj^jjQ ) 
versus halo mass. Our simulated halos are presented (sym- 
bols) alongside the same quantity as it is estimated from 
comparisons of galaxy and halo clus tering and abundances 
dMoster et al.ll2010t IWake et al.1 1201 U curves). The baryon 
conversion efficiency is a strong function of mass in our sim- 
ulations, though not as strong as observations indicate, such 
that small dark matter halos i n our simulations p robably host 
too massive galaxies (see also lDave et al.ll201 lbl) . 

These observed correlations are important constraints that 
should be reproduced in realistic simulated galaxies. The SFR 
versus stellar mass 'main sequence' relation has received sig- 
nificant attention in recent years and it is a key constraint for 
galaxy evolution. The gas fraction of galaxies is a dominant 
factor in their dynamics, and the baryonic conversion effi- 
ciency controls the number density of galaxies. Regardless 
of the imperfect agreement with observations (in particular at 
halo/galaxy masses below the scale that is the focus of this 
paper), the global correlations in Figure Q] are a result of the 
galactic winds model we are using, withou t which they would 
compare much wo rse to the observations (Dave 2008, l2009t 
IDave et aT1l201 lafb|). For example , the three clumpy disks 
simulated by ICeverino et al.l (1201 Ol) in a cosmological con- 
text with less effective feedback have gas fractions as low as 
18%, 14%, 4% and baryon conversion efficiencies as high as 
0.26,0.22, 1.14, both indicating too little suppression of star- 
formation when compared to observations. 

It is apparent that our prescriptio ns dump more momen - 
tum per baryon than used in, e.g., ICeverino et al.l d2010l) . 
iBournaud et al.1 (120071) and lAgertz et all (120091) . As such, 
our sub-grid models may seem extreme and t he ne gative im- 
pact of winds on the clump survival (Section 13.3b overdone. 
However, consider what would happen if we switch off this 
recipe, or substantially turn down its efficacy. To this end, 
we have performed additional experiments, where we used 
less efficient winds, either by injecting less momentum (us- 
ing slower winds, or winds with lower mass-loading fac- 
tors), or by not temporarily decoupling the wind particles 
from the hydrodynamics. If the winds are not decoupled 
they ca nnot escape the vicinity of the galaxies and are inef- 
fective dDalla Vecchia & Schavel2008l) . In these experiments, 
we obtained galaxies with very different properties, in agree- 
ment with similar results from the literature (reviewed above), 
and in contrarst with observations. For example, when halo 
s224 is simulated when winds are not temporarily decou- 
pled from the hydrodynamics, we obtain a central galaxy at 
z = 2 with low specific star formation rate (0.6Gyr~ com- 
pared to 1.5 Gyr -1 in the fiducial model), low gas fraction 
(8.6% compared to the fiducial 40%), high baryonic conver- 
sion efficiency (0.26 compared to the fiducial 0.129), high 
(and peaked) rotation curve (with the peak at 555 kms~ com- 
pared to the fiducial 278kms~ ), and small size (R\/2= 1 kpc 
compared to R\/2 = 3.3kpc). Moreover, the clumpy nature of 
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the galaxies, discussed in detail in the following sections, is 
also lost under such conditions. 

3.2. Disk and clump properties 

In Figure [2] we present gas surface density maps of four 
galaxies from our sample (large panels). The mean gas sur- 
face densities inside radii of pa 4 - 8 kpc have values ranging 
from pa 50-300M Q pc~ 2 . However, the gas surface density is 
far from constant in the disk, but rather exhibits strong local 
maxima, to which we hereafter refer as 'clumps'. The per- 
turbed morphologies, and the clumps, are the result of disk 
instabilities rather than external interactions, as will be dis- 
cussed in Section 13.31 The typical gas surface densities of 
clumps are up to ten times higher than the mean disk sur- 
face densities, i.e. pa 500- 3000 M Q pc -2 . The clumps tend 
to be embedded inside transient ring-like features, although 
not exclusively. The clumpy morphologies of the gas (and 
hence SFR) are clearly visible also after being smoothed 
to the resolution of typical good AO-assisted observations 
(FWHM=1.5kpc), as shown in the small panels, in face-on 
as well as edge-on views (top and bottom small panels for 
each galaxy). After this smoothing, the typical contrast ratio 
(clump-to-disk mean surface densities) is reduced to « 3. The 
total SFR in all clumps is < 20% of the total disk SFR. These 
prop erties match those of observed z sa 2 clumpy disk s very 
well (iGenzel et al.ll201 ltlF5rster Schreiber et al.ll201 lal) . 

We now discuss one galaxy from our sample, s224, in 
greater detail, as a representative case. The gas is distributed 
mainly in a ring with a radius of pa 3 kpc (Figure 0c)). The 
ring, which is transient, formed from spiral features in the 
disk rather than by an interaction (there are no interactions 
with stellar mass ratios < 20 : 1 since z sa 3). Several clumps 
are embedded in the ring, and appear as SFR overdensities 
(Figures |5J a) andQd)). The gas vertical velocity dispersion 
<j z is pa 20- 100 kms -1 (Figure 0b)), in agreement with ob- 
served values at z sa 2. The high density regions and clumps 
are minima of a z , as denser gas dissipates the random motions 
more quickly an d is less prone to s tirring effects such as ex- 
ternal accretion ( lAumer et al.ll2010l) . The disk shows regular 
rotation (Figure [3le )), indicative of its quiet merger history 
(IShapiro et al.ll2008h . 

Hereafter, for s224 we use a threshold of face-on gas sur- 
face density T, g > lOOOM pc~ 2 to define clumps. With this 
definition, we identify four clumps (Figure [2 a)). We calcu- 
late their sizes, masses, gas fractions and SFRs inside the re- 
gion where S ? > 1000 M Q pc -2 , and list them in TableQ] The 
clumps have masses of sa (5 x 10 8 - 1O 9 )M0, U and they are 
very gas rich (/ g ~ 60%). The clumps themselves have low 
(sa 50 kms -1 kpc -1 ) velocity gradients, not much larger than 
the galaxy-wide average velocity gradient, which is typically 
sa 30 kms -1 kpc -1 at a few kiloparsecs from the center. 

1 1 Observed clump masses from the literature tend to be larger, and sev- 
eral clumps with even sa are known. Several biases can contribute 
to this apparent discrepancy. First, our simulated disks are less massive than 
most of the host galaxies of the most massive observed clumps. Therefore, 
it is the clump-to-disk mass ratio that should be compared. Second, some of 
the most massive observed clumps have atypical clump-to-disk mass ratios, 
some of which may be special cases such as clump mergers or galaxy mi- 
nor mergers. Third, there is no standard for identifying clumps and setting 
their borders. Such differences may introduce poorly controlled systemat- 
ics between different studies . The t ypical clump-to- disk mass ratio reported 
in Elmegreen & Elmegreen (2005) is 1%, and in Forster Schreiber et al. 
1201 lbf> it is roughly 2%. This is, at face value, 2 — 3 times higher than the 
clump-to-disk mass ratios found in our simulations. 



3.3. Clump formation and disruption 
Defining a Q parameter for disk instability dToomrdl 19641) 

of 



where a g is the local gas velocity dispersion, c s the gas sound 
speed, and 




(e.g. IDekel et al1 l2009bl) . we find that the clumps are local 
minima with Q <C 1 (Figure [3jf))- The clumps form in-situ: 
they are not accreted from outside the disk, but rather form 
inside the disk via gravitational instability in an environment 
that has Q sa 1 . This is consistent with the theoretical frame- 
work that is described in the Introduction, and in particular in 
agreement with other cosmologic al simulations dAgertz et alj 
l2009tlCeverino et al.l201dl201lb . The clump masses are 1—2 
orders of magnitude larger than the thermal Jeans mass, but 
match the 'largest unstable Jeans scale' (the 'Toomre mass') 
in a turbulent disk 

Mj, a « (^) 2 M disk > 

(|) 2 x6x 10 10 M Q sa4x 1O 8 M (5) 

( IGenzel et all 120081: IDekel et alj l2009bl) . which is driven to 
high values by the large random motions in the disk. 

A notable characteristic of the clumps is their short life- 
times, sa 50Myr, or about half a disk orbital time. The top 
row in Figure @] demonstrates the disruption of a clump in 
s224 by showing a time series of gas surface density maps. In 
the upper right-most panel of Figure[4] the gas masses of three 
clumps are shown as a function of time (solid), demonstrating 
their rapid formation and disruption. In contrast, the bottom 
row of Figure [4] shows the evolution of the same clump in 
an experiment where we temporarily turn the wind off. The 
clump does not disrupt, instead it collapses further and viri- 
alises, and subsequently migrates to the galaxy center. This 
demonstrates that the clumps are destroyed by the wind feed- 
back, and do survive (our equation-of-state is sufficiently soft 
to allow for this), as found in previous simulations, when such 
feedback is absent (or much weaker than adopted in our sub- 
grid model). 

The reason for the disruption of clumps by the wind is 
the following. In our model, the wind mass-loading factor 
(i.e. the outflow r ate over the SFR) for 10 105 M W galaxies 
at z = 2 is r\ sa 4 dOppenheimer & Pavel 120081) . and the ve- 
locity ranges between sa 400 -700 kms -1 . These wind pa- 
rameters agree well with recent obs ervational detection s of 
gas outfl owing from z ~ 2 galaxies dSteidel et alJl20l6T) and 
clumps dGenzel et alj 1201 lb . Star-formation in the clumps 
proceeds o n a ti mescale T$f sa 300 Myr (according to the 
iKennicuti d!998l) relation), thus winds drive gas out of the 
clumps on a timescale sa Tsf/v ~ 100 Myr, which is compa- 
rable to the disk orbital time. The consequence is that the gas 
surface density in the clumps decreases faster than the rate at 
which it is replenished by the instability inside the disk, and 
so clump regions move from Q < 1 to Q > 1 (and thus stop 
collapsing) within a time that is shorter than the disk orbital 
time. The short lifetimes we find are consistent with the upper 
limits of sa 300 -500 Myr derived for the stellar populations 
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TABLE 1 
Galaxy and clump properties 





halo mass 


baryon conversion 


peak 


Half-mass 


Stellar mass 


Gas mass 


Gas fraction 


SFR (b ' 


Number of 


— SFR" Id) 

(SFR) 




M 200 


efficiency 


Vc = ^/GM/R 


radius Rih^"' 


M„< b > 


M g < b > 


f (b) 

Jg 




clumps (c) 




[10 I2 M Q ] 


», 

_()fl b / in 


[kms -1 ] 


[kpc] 


[1O 1O M ] 


[10 10 M o ] 




[Meyr- 1 






median SINS disk M 






237 


4.6 


4.4 


1.9<« 


30% (t) 


65 


5.5 




median simulated disk 


1.5 


0.099 


238 


3.5 


2.4 


2.1 


43% 


27 


5 


1.31 


s224 


1.90 


0.129 


278 


3.3 


4.1 


2.7 


40% 


62 


5 


1.45 


s263 


2.18 


0.097 


289 


3.5 


3.6 


2.1 


37% 


38 


3 


1.51 


s361 


1.29 


0.111 


238 


4.6 


2.4 


2.4 


50% 


27 





0.75 


s377 


1.52 


0.084 


252 


3.3 


2.2 


1.6 


43% 


21 


6 


1.00 


s383 


1.73 


0.115 


253 


3.2 


3.4 


2.1 


38% 


29 


6 


0.86 


s396 


1.50 


0.099 


225 


3.8 


2.5 


1.2 


33% 


14 


5 


0.65 


s447 


1.11 


0.089 


229 


3.5 


1.7 


1.8 


52% 


23 


4 


1.60 


s466 


1.39 


0.094 


238 


5.7 


2.2 


3.0 


58% 


30 


7 


1.41 


s809 


0.83 


0.121 


207 


2.6 


1.7 


1.3 


43% 


19 


5 


1.31 


clump-1 








0.45 


0.022 


0.044 


67% 


2.7 






clump-2 








0.27 


0.02 


0.04 


67% 


2.3 






clump-3 








0.22 


0.022 


0.025 


53% 


1.5 






clump-4 








0.42 


0.027 


0.07 


72% 


3.7 







NOTE. — (a) For galaxies: stellar (3D) half-mass radius within 0.2R200, for clumps: gas half-mass radius inside the region with S g > 1000 Mq pc~ 2 ; (b) 
For galaxies: inside 3R]/2, for clumps: inside the region with S ? > 1000 Mq pc~ 2 ; (c) Identified as local overdensities in the gas surface density, excluding 
the galaxy center ; (d) (SFR) is the mean SFR between z = 3 and z = 2 ; (e) Calculated using the nine galaxies from Forster Schreiber et al. (2009) that are 
defined as 'disks' and have kinematic modeling. Differences to the definitions adopted for the simulations are: V c is the observed maximum rotation velocity 
of the Ha-emitting gas; R\f2 is calculated from projected Hct light; all quantities are galaxy-integrated values; the clumps for MD41, BX389 and BX610 are 
overdensi ties of rest-frame op tical light. Stellar masses and SFRs are esti mated with SEP modeling, and gas masse s are estimated from the SFRs by usin g the 
empirical Bouche et al. (2007) star-formation relation. Data are taken from Genzel et al. 12008); Cresci et al. (2009); Forster Schreiber et al. (2009, 201 la), and 
are all scaled for a[Chabrier ( 2003) init ial mass function ; If) The ga s masses are proba bl y lower bound s , since we use the steep Bouche et al. ( 2007) star-formation 
relation to derive them, as in IForster Schreiber et al. (2009). If the Kennicutt ( 1998) or Genzel et al. 12010) relatio ns were used, t he gas ma sses would increase 
by up to a factor 3. That would increase the gas fractions to 50% — 60%. Indeed, the direct CO observations in Tacconi et al. 12010) and Daddi et al] j2010t) 
indicate 45 — 65% as the most favourable values. 



in the clumps dElmegreen et al.ll2009l: IForster Schreiber et al.l 
1201 lab . This is because the mean stellar ages of the clumps 
in our simulations are larger than the lifetimes of the clumps 
as star-forming overdensities (see FigureQand accompanying 
discussion below). 

Given the masses and sizes of clumps (Table [TJ, their cir- 
cular velocities equal yjGM/R w 80- 100 kms -1 . This is sig- 
nificantly larger than the velocity dispersion of the gas within 
them or their rotational velocities (in particular if the galaxy- 
wide average velocity gradient is subtracted as 'background'). 
Typically the clumps have v 2 0t +2a 2 g rj QMGM/R), or | § | > 3, 
where U is their gravitational potential energy and K the inter- 
nal kinetic energy. Thus, they are neither supported by pres- 
sure nor by rotation, and in fact they are not virialised. This 
is due to continuous feedback on a timescale shorter than the 
dynamical time, preventing virialisation before disruption. 

The clumps consist two types of stellar populations. First, 
'background stars', which are part of the stellar disk at the 
locations where clumps form but are not necessarily bound to 
them, with an age distribution similar to that of the stellar disk 
as a whole. A dynamically cold stellar component also takes 
part in the instability, in addition to the gas, such that the local 
density of old 'background stars' is increased in some clumps 
compared to the clump surroundings. The 'background stars' 
comprise most of the stellar mass of the clumps, as the clumps 
typically convert just about 10% of their peak gas mass into 
stars throughout their lifetime. Second are 'clump stars', 
which form over the lifetime of a clump, and are almost all 
bound to their clumps, constituting typically up to w 20% of 
their stellar masses. 

As a result, the mean age of the stellar populations inside 
clumps follows the global spatial age distribution of the 'back- 



ground stars' of the disk, with an additional small contribution 
of the 'clump stars', which makes the clumps typically be lo- 
cal age minima. We find that the stellar disks in our simula- 
tions have a negative age (and metallicity) gradient, which 
will be discussed in detail in future work. Therefore, the 
mean ages of clumps also present a spatial gradient, such that 
clumps with older stellar populations are closer to the galaxy 
center. This is shown in Figure[5] The trend we find is qualita- 
tively consistent with observations, in spite of the fact that no 
individual clump actually migrates to the galaxy center as its 
stellar populations become older. A quantitative comparison 
to observations cannot be performed yet, as measuring clump 
ages requires high-resolu tion, deep data taken wi t h variou s in- 
struments. For example. IForster Schreiber et alj (1201 lbl) was 
able to measure clump ages for only one z ~ 2 galaxy from 
their sample, BX482, which are shown in Figure|5](j'^Mflre5). 
These measurements rely on Ha equivalent width, and they 
are very sensitive to the assumptions made regarding the star- 
formation history. Therefore, they should be taken only as rel - 
ative ages, not absolute ones (IForster Schreiber et al.ll20lTbh . 
A linear fit in the log-log plane to the clump age-distance re- 
lation from our simulations yields a slope of -0.57 ±0.14. 
The corresponding slope for the BX482 points shown in Fig- 
ure[5]is -2.06 ±0.63. Future comparisons of this correlation 
with larger galaxy samples, and with an improved treatment 
of possible biases both on the theoretical and observational 
sides, may be able to distinguish the s lope w e find from the 
-1.05 slope found by ICeverino et al.l d201 ll) for the case of 
migrating clumps. To conclude, the mere existence of an age 
gradient should not necessarily be in terpreted as supporting 
the clum p migration scenario (see also lForster Schreiber et alj 
d201 lbl) for additional alternative interpretations). 
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FIG. 2 . — Gas surface density maps (in Mq pc~- on a log scale) of four simulated galaxies from our sample. The large panels are maps with resolution compared 
to the spatial resolution of the simulation. The small panels show the corresponding maps after they have been smoothed with two-dimentional Gaussians with 
FWHM=1.5 kpc, in order to imitate the appearance of the galaxies under the highest available resolution of AO-assisted observations. The small panels show both 
face-on (top) and edge-on (bottom) views. All panels are 20 kpc on a side, and all the mag enta bars are 4 kpc long. It is evident that the qualitative apperance o f 
the simulated galaxies is similar to the clump-clusters and chain galaxies observed at z ~ 2 (Elmegreen & Elmegreen 2005, 2006; Forster Schreiber et al. 201 lb). 



4. DIRECT COMPARISON TO SINFONI OBSERVATIONS 

Figure [6] shows mock SINFONI maps for s224 of Ha in- 
tensity, velocity and velocity dispersion. They were obtained 
by converting SFR to Ha l uminosity using Lnatergs -1 ] = 
1 .26 x 10 41 x SFR[ M yr" 1 ] dKennicuttl 1998b . 'placing' s224 
at z = 2.2, convolving it with a Gaussian to imitate a resolu- 
tion of FWHM = 0.17", and pasting it into a real SINFONI 
datacube with pixel scale 0.05" that was obtained from the 
observations of a galaxy from the SINS survey, at a redshift 
where no real galaxy exists. This results in realistic resolution 
and noise properties that correspond to a representative total 
integratio n time (~ 6 hours) for our high-resolution SINFO NI 
data sets dForster Schreiber et al.1120091 iGenzel et al.ll201 ll) . 

The dumpiness, smooth velocity field and relatively flat 
velocity dispersion map (outside the center), which are the 
characteristics of real SINS clumpy disks, are all well repro- 
duced when 'observing' our simulations. The most signifi- 
cant difference to the non-degraded images in Figure [3] (the 
left and middle columns can be directly compared in Figures 
[6]and[3]l is in the velocity dispersion map. The 'beam smear- 
ing' increases the apparent velocity dispersion where there are 
velocity gradients (this is the reason for the diagonal feature 
of high dispersion in the inclined image in Figure |6jc)), and 
the clumps are no longer seen as clear minima. As a result, 
the small but present velocity gradients across the clumps get 



smoothed out. Thus, because of the non-virialised state of the 
clumps, and their being minima in velocity dispersion, they 
do not show strong features in observations at the currently 
available resolution of k 1-2 kpc, even if their masses are 
dynamically sig nificant. This predict ion is con sistent with the 
observ ations of IGenzel et al.l (1201 ll) (but see ICeverino et al.l 
d2011l) for alternative explanations). 

5. SUMMARY AND DISCUSSION 

In this paper we use 'zoom-in' cosmological hydrodynami- 
cal simulations that include a phenomenological model for ef- 
ficient star-formation feedback in the form of fast momentum- 
driven winds with high mass-loading factors, to investigate 
the formation of star-forming disk galaxies at z ~ 2. Our wind 
model reproduces observed low baryon conversion efficien- 
cies and high gas fractions, and we also find a good kine- 
matical and morphological agreement with observations. In 
addition, we present an important numerical result, namely 
that in the presence of such winds, the giant clumps in z ~ 2 
disks, which form in-situ via gravitational instability, may be 
short-lived and unvirialised. This is because they lose mass 
to galactic outflows at a rate that is higher than the rate of 
their dynamical evolution, such that their collapse is halted in 
less than a disk orbital time after their formation, and before 
they reach virial equilibrium. We show that this scenario is 
plausibly consistent with most of the available observations 
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FIG. 3. — Gas properties of the central galaxy s224 at z. = 2. (a): face-on SFR surface density, (b): face-on line-of-sight (vertical) velocity dispersion, (c): 
face-on gas surface density, (d): edge-on SFR surface density, (e): edge-on line-of-sight velocity field overplotted with contours of SFR surface density and (/): 
log Q. Each panel is 10 kpc on a side. All face-on panels are overplotted with E g = 1000 Mq pc~ 2 contours. 




FIG. 4. — A time sequence of gas surface density maps showing the disruption of a clump in our model (Top), where t = (not shown) is the formation time 
of the clump. To demonstrate the role of the wind, we turn it off at z = 2.03 (t = 22 Myr) and show the alternative evolution of non-disruption, virialisation and 
migration (Bottom). The color coding is as in Figure 0c). The upper right-most panel shows the mass of gas (solid) and young (< 50 Myr) stars (dashed) for 
four clumps as a function of time since their formation. The magenta lines are for the clump highlighted on the Top and the black for the clump highlighted on 
the Bottom. The jump in mass of the green lines at t « 60 Myr is a result of a merger between two clumps (not shown in other panels). The typical clump lifetime 
in the presence of winds is es 50 Myr, and the mass of new-formed stars is approximately 10% of the maximum clump gas mass. The mass of new-formed stars 
internal to the clump decreases following the decrease of the gas mass, as these stars are dispersed out of the clump when the collapse is halted by the return to 
Q> 1. 
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FIG. 5. — Clump stellar age (mass-weighted) versus distance from the 
galaxy center. Colors correspond to the different galaxies in the same way as 
in Figure [T] both for z = 2 (dots) and z = 2.2 (circles). More distant clumps 
tend to be younger, because of the dominance of 'background' stars in the 
clumps, and given the global age gradient of the galaxy. A similar relation 
for the clumps in the observed galaxy BX482 ('raw' measurements, without 
backgr ound subtraction) is shown as well (squares; [Fbrster Schre iber et all 
I201 lbf) . Given that the observed ages are only relative ages, our simulations 
seem consistent with the observations of this galaxy. The important conclu- 
sion is that the observed trend of clump stellar age with distance (see also 
Genzel et al. 2008; Elmeareen et al. 2009) is not necessarily an indication for 
clump migration to the center of the galaxy, as has been previously consid- 
ered. 



of clumpy disks at z ~ 2: their global properties (Figure [TJ, 
their dumpiness (Figures [2] and [3), their kinematics and ap- 
pearance in SINFONI observations (Figures[3]and[6]i, and the 
gradient they show between clump stellar age and clump dis- 
tance from the galaxy center (Figure |5}. 

Our results are in contrast to previous models that do not in- 
clude fast winds with high mass-loading factors. We demon- 
strate that if we temporarily turn off the winds, we reproduce 
previous results where clumps are long-lived and migrate to 
the galaxy centers. Therefore we are able to directly show that 
the winds are the critical component of our model that makes 
the clumps short-lived. However, if the wind is drastically re- 
duced starting from the initial conditions of the simulations 
at high redshift, gas-rich unstable disks do not develop in the 
first place. Instead, the galaxies become gas-poor, very mas- 
sive (compared to the halos inside which they form, i.e. with 



very high baryon conversion efficie ncies), and they d evelop 
very peaked rotation curves (see also lJoung et al.ll2~009h . 

It should be stressed that our model of the winds is phe- 
nomenological, and our results are based on the working hy- 
pothesis of the existence of powerful superwinds with some 
given, a priori properties. Our hypothesis still awaits fur- 
ther verification. For example, our feedback recipe assumes 
that gas can be blown out of the clumps and subsequently 
out of the galaxy at an imposed rate proportional to the star- 
formation rate. It remains to be better understood how this 
occurs on sub-clump scales b ased on improved physical mod- 
els (e.g. lHopkins et al1l201 ll) . Also, the momentum we inject 
into the wind per unit stellar mass formed is at the high side 
of the theoretical expectations. While winds are observed to 
be ubiquitous in high-redshift galaxies, and in particular they 
have been recently observed from individual giant clumps, 
the measurements of their crucial parameters, in particular the 
mass-loading factors, is still very uncertain. 

It has been suggested in the recent literature that the migra- 
tion of giant clumps to the centers of high-redshift galaxies 
may serve as an avenue for the formation of galaxy bulges, 
in a ddition to the more tradition al avenue of galaxy mergers 
(e.g. lSteinmetz & Nav arro 2002]). The alternative scenario we 
present in this paper, namely that giant clumps do not migrate 
to their galaxy centers, does not necessarily mean that bulges 
do not form secularly in high-redshift galaxy disks. The in- 
stabilities and irregularities we find in our simulated disks do 
induce internal torques and angular momentum loss, such that 
inflows to smaller radii inside the disks do occur. Neverthe- 
less, by comparing clump evolution with and without winds, 
we demonstrate that secular bulge growth occurs much more 
slowly in the case where winds disrupt clumps. We reserve 
detailed studies of secular bulge formation under conditions 
of strong feedback and giant clump destruction, as well as 
implications for the global cosmic rate of bulge formation, to 
future work. 
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